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Abstract 


This  Report  expands  upon  our  previous  work  in  this  area.  Specific  appli¬ 
cation  is  made  to  searches  for  near-stationary  artificial  satellites  and  to 
the  classical  along  orbit  search.  These  two  reflect  different  limiting  cases 
for  the  a  priori  target  distribution  (uniform  for  the  near-stationary  case) 
and  the  conditional  detection  probability  (uniform  for  the  along  orbit  case). 

Our  treatment  of  the  near-stationary  case  is  as  realistic  as  Is  currently 
possible.  Atmospheric  extinction,  eclipses,  and  phase  effects  are  all  Included. 
Similarly  we  have  explored  a  variety  of  scenarios  for  the  along  orbit  search. 

We  conclude  with  explicit  search  plan  construction  and  illustrative  examples. 
Finally  we  mention  other  work  in  this  area  currently  underway. 
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I.  REVIEW 


In  reference  I  one  of  us  (LGT)  introduced  the  concepts  of  optimal  search 
theory  to  the  artificial  satellite  search  problem.  The  best  summary  of  the 
literature  (reference  2)  is  herein  extended  in  detail  to  two  classical  search 
problems  fo1"  deep  space  artificial  satellites--near-stationary  satellites  and 
along  orbit  searches.  Clearly  the  latter  type  Is  frequently  utilized  in  the 
low  altitude  regime  and  by  radars  as  well  as  by  passive  optical  sensors. 

Hence  the  formalism  that  allows  one  to  solve  such  problems  should  be  of 
interest  to  a  wide  audience. 

A.  Formulation 

One  looks  for  artificial  satellites  on  the  celestial  sphere.  In 
the  largest  sense  this  forms  the  two  dimensional  search  space  of  the  problem. 

In  practice  we  delineate  a  limited  area  of  the  celestial  sphere  (say  above 
altitude  30°)  that  we  shall  actually  search  In.  Denote  this  search  space  by  J. 

One  searches  using  a  telescope  with  a  finite  field  of  view.  In  practice 
we  always  examine  an  entire  field,  never  a  fraction  of  a  field  nor  more 
than  one  field  at  a  time.  Hence  the  search  space  J  consists  of  a  discrete 
set  of  fields  of  view.  Number  these  by  the  index  j  =  1,  2,  .  .  .  In 
particular,  since  the  celestial  sphere  encompasses  4tt  steradians,  max  (j)  <  ®. 

Before  the  satellite  is  found  one  assigns  an  a  priori  target  distri¬ 
bution  on  the  search  space  J,  p:  J  -*■  [0,1]  (the  notation  means  that  p  Is  a 
function  defined  over  the  set  J  which  maps  elements  of  J  into  the  domain 
zero  to  unity  Inclusive).  The  target  distribution  is  the  a  priori  proba¬ 
bility  of  finding  an  artificial  satellite  in  field  of  view  j  e  J  before  one 
starts  the  search.  For  near-stationary  satellites  a  reasonable  model  for  p 
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Is  p  is  uniform  over  all  geocentric  right  ascensions  and  over  the  geocentric 
declination  range  <  10°  (or  5°  or  20°).  For  other  satellites,  both  because 
of  parallax  effects  and  the  Inherent  spread  over  orbital  element  space 
(particularly  In  Inclination),  a  reasonable  model  for  p  is  that  p  is  uniform 
over  the  topocentric.  celestial  sphere.  In  any  case 

£  p(j)  <  1 
JeJ 

When  one  examines  a  field  of  view  for  an  artificial  satellite  one  expends 
a  certain  amount  of  effort  trying  to  detect  the  satellite.  One  may  look  at 
the  same  field  of  view  several  times.  The  cost  of  performing  k  inspections 
In  the  j'th  field  of  view  is  measured  by  a  cost  function 

c(j,k):  J  x  {0,1,2,...}  -*•  [OH.  Clearly  c(j,0)  =  0  V  j  e  J  (no  effort  implies 
no  cost).  One  could  measure  cost  by  the  time  spent  examining  a  field  of  view 
plus  the  time  spent  in  moving  to  the  next  field  of  view  (this  makes  c  non¬ 
local  and  is  not  desirable).  Operationally  we  always  spend  the  same  time  in 

1  /2 

each  field  of  view  (mere  or  less).  Also,  because  [area  (0)]  /slew  speed  « 
time  spent  examining  a  field  of  view,  the  non-local  element  of  c  is  both 
unimportant  and  varies  little.  Thus  we  shall  measure  cost  by  time  and 
specialize  to  the  case  when  the  Incremental  cost  of  the  k'th  examination  in 
field  of  view  number  j,  viz. 

y(j,k)  -  c(j,k)  -  c(j,k-l) 

is  a  constant  Independent  of  both  j  ( i . e . ,  the  telescope  is  fast  and  all 
fields  of  view  are  treated  equally)  and  k  (e.g.,  the  same  field  of  view  is 
equally  well  Inspected  each  time). 


When  one  does  examine  a  field  of  view  of  the  search  space  looking  for 
an  artificial  satellite  then  there  is  a  conditional  probability  of  detecting 
it  on  or  before  the  k'th  inspection  of  that  field  of  view  (given  that  it  is 
there).  This  function,  for  field  of  view  number  j  and  examination  k,  is 
denoted  by  b{j,k):  J  x  {0,1,2,...}  -*•  [0,1].  Naturally  b(j,0)  =  0  V  j  e  0 
(you  can't  find  it  if  you  don't  look  for  it).  From  the  detection  function  b 
one  can  construct  the  probability  of  failing  to  detect  the  satellite  on  the 
first  k-1  scrutinizations  of  field  of  view  number  j  and  then  succeeding  on  the 
k'th  one  (given  that  the  satellite  is  in  field  of  view  number  j);  viz. 

3(j ,k)  =  b(j,k)  -  b(j,k-l) 

There  is  a  lot  of  physics  and  mathematics  subsumed  in  the  detection 
function.  Clearly  It  depends  on  the  satellite's  apparent  magnitude,  the 
background  star  density,  the  night  sky  background  brightness,  the  resolution 
element  size  of  the  detector(s),  the  false  alarm  probability  one  is  willing  to 
accept,  how  tired  one  is,  etc.  Since  the  celestial  sphere  is  unchanging, 
atmospheric  extinction  can  be  computed,  the  Moon's  position  is  known,  etc. 
this  Is  a  computable  function.  Operationally,  for  a  fixed  set  of  external 
parameters,  our  detection  probability  has  the  shape  shown  in  Fig.  1  where 
m^  is  our  quoted  limiting  magnitude  (e.g.,  where  the  probability  of  detection 
Is  50J).  The  form  shown  In  the  diagram  will  be  used  to  compute  the  optimal 
search  plans  given  below. 

Finally  we  need  to  define  a  search  plan.  A  discrete  plan  is  a  sequence 
K  -  (£]  *  *  * )  wblch  tells  the  searcher  to  first  look  in  cell  C-j ;  if  the 

satellite  is  found  there  then  terminate  the  search  but  if  it  isn't  found  then 
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look  next  in  field  of  view  etc.  A  global  way  to  describe  this  is  by  a 
function  which  specifies  the  allocation  of  effort  devoted  to  each  field  of 
view  j.  To  this  end  define  f:  J  -*•  [0,®);  f(j)  is  the  number  of  examinations 
in  field  of  view  number  j. 

Above  we  referred  to  searches  for  a  fixed  target.  Clearly  the  satellites 
we  are  trying  to  find  are  moving.  We've  made  the  assumption  that  these  objects 
are  fixed  when  compared  to  our  search  rate.  The  mathematical  formulation  of 
this  approximation  is  [area  (J)/search  rate] -satel 1 1 te  angular  speed  «  field 
of  view. 

B.  Optimal  Searches 

Given  the  cost  of  searching  field  of  view  number  j  a  total  of  k 
times,  c(j,k),  the  total  cost  of  performing  the  search  plan  £  with  allo¬ 
cation  f  is 

C[f]  *  c(j,f(j)) 

The  total  number  of  examinations  over  all  fields  of  view  is  ~ 

jeJ 

Similarly  the  total  probability  of  satellite  detection  with  this  allocation  of 
effort  is  P[f], 

P[f]  -  £  p( j )b( j ,f( j ) ) 
jeJ 

where  b(j,k)  Is  the  conditional  probability  of  finding  the  satellite  In  field 
of  view  number  j  after  k  examinations  of  that  field  of  view  given  that  it’s 
in  that  field  of  view. 
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There  are  four  types  of  searches  one  might  define  as  optimal.  One  might 
be  interested  in  maximizing  the  total  probability  of  detection  when  constrained 
to  a  given  number  of  inspections  {say  K).  If  the  incremental  cost  function 
y(j.k)  =  c(j,k)  -  c(j,k-1)  is  a  constant,  then  (after  a  suitable  renormaliza¬ 
tion)  one  is  demanding  that  P[f]  be  a  maximum  for  C[f]  _<  K.  Such  a  search  is 
termed  totally  optimal .  If  one  demanded  optimality  for  all  K  =  1,  2,  3,  .  .  . 
then  the  search  is  called  uniformly  optimal .  A  third  type  of  search  plan 
that  one  might  consider  is  the  search  plan  that  maximizes  the  probability  of 
detection  with  respect  to  the  incremental  cost  and  does  so  at  every  step  of 
the  search.  Mathematically  one  finds  the  value  of  j  which  maximizes 
p(j)B(j,k)/Y(j»k)  at  each  k.  These  searches  are  called  locally  optimal . 

Lastly  one  might  entertain  a  search  plan  that  minimized  the  total  expected 
cost  ( 1 . e . ,  was  the  fastest)  to  find  the  target. 

The  essential  assumptions  necessary  to  cast  the  artificial  satellite  search 
into  the  simplest  form  of  the  mathematical  superstructure  that  Stone2  outlines 
are 

(1)  That  the  satellite  is  fixed  (l.e.,  search  rate  high 
compared  to  the  satellite's  angular  speed), 

(2)  That  the  search  space  is  discrete  (i.e.,  a  fixed  field  of 
view) , 

(3)  That  the  allocation  of  effort  is  discrete  (i.e.,  no  favored 
fields  of  view) ,  and 

(4)  That  y  Is  bounded  away  from  zero  and  p( j )b( j ,k)/Y( j  ,k)  is  a 
decreasing  function  of  j  (i.e.,  no  free  examinations  of  a 

field  of  view  and  the  larger  the  search  space  the  more  difficult 
to  detect.) 
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I  do  not  believe  that  the  physics  or  astronomy  is  strained  by  these  strictures 
In  fact  (5)  y  =  constant  is  not  unreasonable  (i.e.,  the  telescope  moves  smart¬ 
ly).  The  important  point  is  that  under  these  five  limitations  the  totally 
optimal  search  plan,  the  uniformly  optimal  search  plan,  the  locally  optimal 
search  plan,  and  the  fastest  searches  are  all  identical .  Not  only  that,  it 
can  be  explicitly  exhibited.  See  Stone's  text  for  the  rigorous  mathematical 
statements  of  the  relevant  theorems  and  their  proofs. 

C.  The  Search  Plan 

We  need  just  a  bit  more  mathematics  before  we  can  exhibit  the  solu¬ 
tion  to  the  optimal  search  problem.  The  search  plan  £  =  (£, ,  £2»  £3*  •••)  i<: 
a  sequence  of  values  Cj  e  J  for  i  =  1,2,3,  ...  These  specify  th^t  the  i'th 
examination  be  in  field  of  view  E,.  if  the  previous  i-1  inspections  failed  to 
detect  the  satellite  in  fields  of  view  £^,£2,...,  Let  the  set  of  all 

such  search  plans  be  denoted  by  E.  Introduce  the  probability  P[n,  £]  (and 
the  cost  C[n,  £])  of  detecting  the  satellite  on  or  before  the  n'th  ex-mi- 
nation  while  performing  search  plan  £  c  E  (of  the  first  n  inspections). 
Finally,  let  r(j,n,£)  be  the  number  of  scrutinizations  out  of  the  first  n 
that  are  placed  in  j'th  field  of  view  while  following  search  plan  £.  A  uni¬ 
formly  optimal  search  plan  [for  y(j,k)  =  1;  this  is  an  unimportant  normaliza¬ 
tion]  £*  e  E  is  one  such  that 

P[n,£*]  =  max  {P[n,£]  :  £  e  E  },  n  *  1,  2,  .  .  .,  K 

A  locally  optimal  search  plan  £*  is  one  such  that  £^  is  determined  by 
[y  f  0  necessarily] 


P(C1)B(?1,1) 
y(«,,D  j£j  r(j,D 

and  having  determined  the  field  of  view  for  the  first  n-1  examinations 
(£.j,  ])  field  of  view  for  the  n'th  one  is  determined  from 

p(i)B(i,r(i,n-i,0  +D  max  p(j)B(j,r(j,n-1  ,£)  +1) 
y(i ,r(i ,n-l ,C)  +  1)  JeJ  y(j ,r( j ,n-1 ,£)  +  1) 

with  Cn  =  i.  Now  define  kn  =  r(Cn>  n,  £)•  The  notation  means  that  the  n'th 
examination  of  the  search  plan  5  is  placed  in  field  of  view  £n  and  that  it  is 
the  kn'th  time  that  this  field  of  view  has  been  searched.  The  average  cost  to 
find  the  satellite  can  be  expressed  in  a  variety  of  ways  if  the  limit  as 
n  -►  oo  of  P[n,5]  is  unity; 

CO 

y(C)  =  E  C[n,d  (P[n,5]  -  P[n-l,C]) 

n=l 

«°  n 

■  Y<Vkn>  P<«n>  6<Vkn> 

n=l  m= l 

“  £  £  Y^m’kn^ 

m= l  n-m 

CO 

=  E  Y(^fn.km)  (1-P[m  ,cj) 

m=l 

since  P[0,5]  =  0.  If  Y(j»k)  =  1  then  this  reduces  to 

y(0  =  E  ( 1  -P[n  ,C] ) 
n=0 
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Now  we  can  exhibit  the  solution  explicitly.  Under  the  assumptions  out¬ 


lined  above  if  qj  Is  the  probability  of  detecting  the  satellite  after  a  single 
examination  of  field  of  view  number  j  {given  that  it  is  in  field  of  view 
number  j)  then,  as  each  inspection  is  an  independent  event,  the  incremental 
conditional  probability  of  detection  8 ( j  *  k )  =  b(j,k)  -  b(j,k-l)  is  given  by 

g( j ,k)  =  qj  (l-qj)k_1  for  j  e  J  ,  k  =  1,  2,  .  .  . 

Normalize  such  that  y(j,k)  -  1  v  j  e  J,  k  =  1 ,  2,  .  .  .  and  suppose  that  an 
allocation  f(j)  has  total  cost  (i.e.,  number  of  inspections)  K, 

£  f(j)  =  K 
jeJ 

The  total  probability  of  detection  for  this  allocation  of  effort  will  be 

~  v,  f(j) 

P[f]  =  L  p(j )b(j ,f ( j ) )  =  E  p(.i)[l-(l-qi)  ] 

jcJ  jeJ  J 


Consider  the  search  plan  defined  by:  one  makes  the  n'th  inspection  1r 
field  of  view  number  1  e  J  such  that 


r(i ,  n,  O  r(j,  n-1 ,  £) 

p( i )q1  ( l -q1 )  =  p(j)qj  U  -  qj) 

Then  C  *  C  and  is  optimal  (in  all  senses).  This  result  is  due  to  Chew^. 
Since  J  is  finite  the  existence  of  an  i  satisfying  the  above  is  guaranteed. 
If  one  exploits  the  uniformity  of  the  target  distribution  p  over  the  search 
space  J,  then  the  result  is  even  simpler. 


-  q4) 


r(i,  n,  O 


max 

jeJ 


q  j  ( 1  51 


0) 
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II.  NEAR-STATIONARY  SATELLITE  SEARCHES 

We  have  already  argued  that  the  a  priori  target  distribution  p(j)  can 
be  approximated  by  a  defective  uniform  distribution  over  the  search  space. 

We  have  also  argued  that  the  Incremental  cost  function  is  homogeneous  over  J 
and  Independent  of  the  numbe"  of  looks,  y(J*k)  •  1  (In  appropriate  units). 

The  probability  of  detection  of  the  satellite  In  field  of  view  number  J  (given 
that  It's  there)  Is  dj.  This  depends  principally  on  the  apparent  magnitude  of 
the  satellite  and  the  night  sky  background.  Three  effects  tend  to  make  satel¬ 
lites  fainter;  atmospheric  extinction,  loss  of  brightness  due  to  Increasing 
phase  angle,  and  Increasing* distance  (heliocentric  or  geocentric). 

The  extinction  Is  modeled  a*  usual, 

t  "  cz  sec z 

where  z  Is  the  topocentrlc  zenith  distance  and  c2  Is  the  extinction  per 
unit  air  mass.  We've  used  a  value  of  0.13  mag/air  mass  for  cz>  For  the  phase 
function  In  magnitudes  I've  used  the  results  In  reference  4 

B(1 ,0)  -  B(1,0)  ♦  0.530  -  0.134  |0|0*714  -  71.  for  |0|  <  7° 

B(1 ,0)  •  8(1 ,0)  -  et  for  | 0 |  2  7° 

where  8(1,0)  Is  the  absolute  B  magnitude  and  B(l,9)  Is  the  apparent  magnitude 
corrected  for  phase  angle  0.  The  parameter  of  the  linear  part  of  the  phase 
function  In  magnitudes  l  ■  0.039  mag/deg. 


Eclipses  can  play  an  Important  role  in  searches  for  near-stationary 
satellites.  Obviously  it  makes  no  sense  to  search  that  part  of  the  sky  cur¬ 
rently  undergoing  eclipse  (in  the  visible  wavelength  bands).  Moreover,  it 
would  be  a  mistake  to  preferentially  find  an  artificial  satellite  just  prior 
to  its  eclipse  (eg.  westward  of  the  shadow  cone).  Hence  not  looking  within 
the  penumbra  represents  one  level  of  search  planning  refinement  and  only  look¬ 
ing  eastward  of  the  penumbra  represents  yet  another  step  up  in  sophistication. 

The  geometrical  discussion  of  artificial  satellite  ecliDses  is  more  com¬ 
plicated  than  that  for  the  Moon  because  the  satellites  are  much  closer.  The 
results  were  worked  out  by  one  of  us  (LGT)  In  1981.  For  a  near-stationary 
satellite  the  half-angle  of  the  umbral  cone  is  8?427  while  that  of  the 
penumbral  cone  Is  8?960.  Since  it  only  takes  ^2  minutes  to  traverse  this 

and  the  brightness  variation  during  the  transition  through  penumbral 
eclipse  Is  exceedingly  difficult  to  model,  we've  chosen  to  ignore  the  distinc¬ 
tion.  The  half-angle  of  the  penumbral  eclipse  cone  was  then  increased  by  2% 
to  allow  for  refraction  effects  In  the  Earth's  atmosphere.  El  1 1 ptl cl ty  of 
the  Earth,  Its  atmosphere,  and  Its  heliocentric  orbit  have  either  been  ignored 
(the  first  two)  or  averaged  over  (the  latter). 

Since  the  eclipse  Is  centered  at  the  solar  opposition  point  we  merely  test 
for  the  angular  distance  from  the  center  of  a  field  of  view  to  the  opposition 
point.  If  this  distance  Is  less  than  1.02  times  the  half-angle  of  the  penum¬ 
bral  shadow  cone  then  the  apparent  magnitude  Is  set  equal  to  ®  which  results  in 
a  probability  of  detection  of  zero  (cf.  Fig.  1).  Otherwise  the  apparent  magni¬ 
tude  Is  calculated  as  Indicated  above.  For  the  trailing  shadow  cone  searches, 
the  right  ascension  of  the  center  of  the  field  of  view  is  less  than  (mod  24) 
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that  of  the  opposition  point,  then  the  probability  of  detection  =  0.  Other¬ 
wise  It  Is  computed  as  just  described.  Sample  search  plans  are  shown  in 
Figs.  2-8. 

All  of  the  search  plans  have  been  terminated  under  the  same  (artificial) 
criteria:  Either  each  field  of  view  of  the  search  space  (30°  x  3h)  In  declin¬ 
ation  x  right  ascension)  Is  examined  once  or  a  field  of  view  of  the  search 
space  Is  about  to  be  examined  for  the  seventh  time  (a  combination  storage/ 
reasonableness  criteria  of  futility).  The  fields  of  view  are  each  2°  x  2° 
(uncorrected  for  the  cos6  foreshortening) .  Furthermore  all  of  the  sample 
search  plans  are  Illustrated  In  the  same  format:  Declination  Increasing  up 
(North)  and  right  ascension  to  the  right  (East).  The  numbers  In  the  Indivi¬ 
dual  fields  of  view  are  the  examination  number,  following  the  optimal  search 
plan,  of  that  field  of  view.  No  entry  means  that  the  field  of  view  wasn't 
examined  before  termination.  Thus,  if  you  look  at  Fig.  2,  the  most  north¬ 
easterly  field  of  view  was  never  examined  while  the  field  of  view  just  to  Its 
west  was  Inspected  tw1ce--on  the  1276  look  of  the  optimal  search  plan  and  on 
the  1372  look.  Similarly  the  highest  probability  field  of  view  was  searched 
on  the  first  examination  of  the  optimal  search  plan  and  on  the  562  look.  The 
effort  Is  distributed  In  elliptical  waves,  about  the  opposition  point.  In 
Fig.  2,  The  major  axis  of  the  ellipse  Is  vertical--jo1n1ng  the  opposition 
point  and  the  zenith.  Actually  the  "ellipses"  are  not  North-South  symmetri¬ 
cal,  they  bulge  more  In  the  North  (eg.  note  the  number  of  empty  fields  of 
view  In  the  southwestern  and  southeastern  corners  as  opposed  to  the  corre¬ 
sponding  northern  ones).  The  reason  for  this  is  that  the  satellites  are 
brighter  nearer  the  zenith  (because  of  reduced  atmospheric  extinction)  than 
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130875-N 


IIBIIHIH 


EAST 


RIGHT  ASCENSION 


Fig.  2.  Optimal  search  plan  for  a  near-stationary  artificial  satellite  on  mid¬ 
night  local  time  of  a  winter  solstice  night.  Each  rectangular  field  of  view  is 
2°  X  2°.  See  the  text  for  a  fuller  description. 

13 


they  are  closer  to  the  nadir.  The  cumulative  probability  of  detection  at  the 
termination  of  the  search  was  60. 7%  assuming  that  the  artificial  satellite  was 
in  the  search  space  initially.  After  examining  J  fields  of  view  the  optimal 
plan  was  23%  more  efficient  than  the  uniform  (eg.  existing)  search  plan.  This 
Increase  in  efficiency  had  dropped  to  4.2%  after  4J  cells  had  been  examined 
(that  Is  the  cumulative  probability  of  detection  of  the  optimal  search  plan 
after  the  firt  4J  examinations  was  4.2%  larger  than  the  cumulative  probability 
of  detection  of  4  repetitions  of  the  plan  that  searches  each  of  the  J  fields 
of  view  of  the  search  space  once  [=  the  uniform  plan]).  This  loss  of  effec¬ 
tiveness  is  easy  to  understand  by  actually  looking  at  Fig.  2.  As  the  search 
wears  on  the  allocation  of  effort  approaches  uniformity,  hence  the  relative 
advantage  should  approach  unity. 

The  next  figure  (Fig.  3)  shows  the  same  scenario  but  at  3  A.M.  local  time 
Instead  of  at  midnight.  There  Is  an  obvious  northeasterly  shift  of  effort 
(because  of  extinction).  There  are  further  small  scale  differences  but  more 
Importantly  after  J  field  of  view  examinations  the  optimal  search  is  34.4% 
more  efficient  than  a  uniform  one  would  be  (with  the  same  total  cost  of 
course).  Figure  4  Is  Identical  to  Fig.  2  but  for  a  satellite  a  half  magnitude 
brighter.  Now  each  cell  of  the  search  is  examined  because  it  makes  sense  to 
do  so.  The  final  cumulative  probability  of  detection  Is  83.1%. 

Figure  5  repeats  the  scenario  of  Fig.  2  but  now  eclipses  are  Included. 

The  hole  In  the  center  represents  the  penumbral  shadow  cone  and  this  search 
plan  Is  markedly  different  from  the  one  In  Fig.  2.  This  time  the  optimal 
search  plan  is  65.7%  more  efficient  than  the  uniform  one  after  J  looks,  42.7% 
after  2J  looks  and  32.0%  after  3  looks.  The  cumulative  probability  of  detec¬ 
tion  at  termination  is  35.9%. 
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Fig.  Same  as  Fig.  2  but  for  a  satellite  of 


mc 
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brighter. 


16 


DECLINATION 


WEST  EAST 

RIGHT  ASCENSION 


Fig.  5*  Same  as  Fig.  2  but  including  eclipses. 
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In  Fig.  6  we’ve  changed  from  the  winter  solstice  to  the  vernal  equinox. 
Otherwise  it’s  a  repetition  of  the  Fig.  3  scenario  but  at  9  P.M.  local  time. 


The  North-South  asynmetry  is  especially  evident  because  opposition  is  lower  in 
the  sky  (6  =  0°)  than  It  was  at  the  winter  solstice  (6  =  e  =  obliquity  of 
the  ecliptic  =  23? 5 ) .  After  an  expenditure  of  effort  of  J  scrutinizations  the 
optimal  search  plan  Is  60.3%  more  efficient  than  is  the  uniform  one.  Figure  7 
is  also  on  riarch  21  but  at  midnight;  compare  with  Fig.  5.  Finally,  Fig.  8 
shows  vernal  equinox,  midnight,  trailing  edge  shadow  cone  search  plan.  The 
first  field  of  view  is  displaced  northeast  (extinction)  and  there's  a  marked 
tendency  to  look  in  the  middle  declination  fields  of  view  (phase  loss). 

Hopefully  these  examples  will  convince  you  that  optimal  search  planning 
is  non-trivial,  non-intuitive,  and  Important. 
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Fig.  6.  Same  as  Fig.  2  but  at  9  P.M.  local  time  on  the  vernal  equinox. 
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Fig.  7.  Same  as  Fig.  5  but  on  the  vernal  equinox. 
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Fig.  8.  Midnight  local  time,  vernal  equinox  optimal  search  plan  including 
eclipses  but  only  trailing  shadow  cone  edge. 
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III.  ALONG  ORBIT  SEARCHES 

One  usually  Initiates  an  along  orbit  search  when  an  examination  of  the 
nominal  position  for  the  artificial  satellite  doesn't  yield  the  object  of 
Interest.  In  addition  one  has  assumed  no  change  in  the  satellite's  orbital 
plane,  merely  that  the  satellite  Is  slightly  late  or  slightly  early.  Hence 
an  along  orbit  search  Is  one  wherein  a  series  of  fields  of  view  (or  beam 
widths),  along  the  topocentrlc  path  of  the  assumed  orbital  plane  and  travel¬ 
ling  with  the  satellite's  presumed  mean  motion,  are  examined  under  the  hypo¬ 
thesis  that  the  mean  anomaly  at  epoch  is  slightly  larger  or  smaller  than  its 
nominal  value.  Such  a  procedure  renders  the  moving  target  fixed  in  the  moving 
coordinate  system  of  the  search  space.  As  the  usual  reason  for  the  non- 
appearance  of  a  satellite  is  a  perturbation  in  the  mean  motion  (atmospheric 
drag)  or  a  poor  elenent  set  (due  to  a  mean  motion/eccentricity  swap),  the 
search  space  is  not  strictly  comoving  with  the  satellite.  Therefore  to  remain 
within  the  fixed  target  scenario,  the  strong  inequality 

area  search  space  <<;  field  of  view 

search  rate  relative  angular  speed 

must  hold  (relative  angular  speed  =  jtrue  mean  motion- nominal  mean  motion)). 

I  shall  assume  this  to  be  the  case. 

With  this  point  In  mind  the  search  snace  J  for  the  along  orbit  search  is 
a  linear  series  of  field  of  view  (see  Fig.  9).  I  shall  label  them  by  j  =  -N, 
-N  ♦  1 ,  ... .,  0,  . ..,  N.  The  a  priori  target  distribution  p  will  be  assumed  to 
be  synmetrlcal  [p(n)  =  p(-n)  n  =  0,  1,  2,  . . . ,  N]  and  unimodal  [l^p(0)>p(  1  )_> 

. ..>p(N)].  We  Intend  to  acquire  data  this  Fall  at  the  ETS  to  provide  an 
empirical  estimate  for  p  as  a  function  of  orbital  type  and  age.  For  the 
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Fig.  9-  Comovlng  and  actual  fields  of  view  for  an  along  orbit  search. 


explicit  numerical  computations  listed  below  we  shall  further  assume  that  the 
probability  of  being  In  either  of  the  next  outer  fields  of  view  Is  equal  to 
that  of  being  In  the  adjacent  Inner  one,  viz. 

P(n  +  1)  +  p(-n  -  1)  =  p(n)  =  p(-n) 

and  since  p  Is  symmetrical  this  Implies  that 

p(n  +  1)  =  p(n)/2  n  =  1,2,  ....  N 

The  cost  function  Is  most  realistically  defined  in  terms  of  t1me--the 
time  spent  examining  a  field  of  view  plus  the  time  required  to  move  to  the 
next  one.  In  the  case  of  the  6E00SS  network,  the  telescopes  accelerate  and 
decelerate  so  fast  that  the  non-local  nature  of  the  cost  function  for  the 
small  areal  extent  (<(2N  +  1)0  where  0  is  the  field  of  view  [N.B.  This  is  not 
strictly  true  because  J  is  comoving  with  the  satellite.  The  real  maximum 
angular  extent  Is  a  complicated  function  of  N,e,  the  satellite's  mean  motion, 
and  the  time  spent  Inspecting  each  field  of  view]}  Involved  here  shall  be 
neglected.  Since  the  time  spent  examining  each  field  of  view  is  a  constant, 

c(j,k>  =  k(tl00k  *  vjd,  k  •  1,2,... 

where  t^  ^  Is  the  time  spent  looking  in  a  field  of  view  for  the  artificial 

satellite  and  t _ _  is  an  average  duration  of  a  movement  from  one  field  of 

move 

view  to  another.  In  the  appropriate  set  of  units  the  incremental  cost  func¬ 
tion  b(j,k)  *  c(j,k)  -  c(j,k  -  1)  =  l. 

The  last  quantity  needed  to  specify  the  problem  is  the  conditional 
detection  probability  b(j,k).  As  the  portion  of  the  celestial  sphere  covered 


by  J  is  small,  the  natural  background  and  foreground  sources  of  noise  vary 
smoothly,  and  the  time  to  complete  an  along  orbit  search  is  short,  the  differ¬ 
ential  extinction,  phase  angle,  night  sky  background,  etc.,  effects  are  all 
small.  Therefore,  the  essential  approximation  for  the  along  orbit  search  is 
that  b(j,k)  is  homogeneous, 

b(i,k)  =  b(j,k)  vi.jcJ;  k  =  1,2,... 

If  b  is  the  single  glimpse  probability  (i.e.,  the  conditional  probability  of 
detection  in  any  field  of  view  given  that  the  satellite  is  in  that  field  of 
view)  then  the  Incremental  probability  of  detection  on  the  k'th  examination 
is 

8(j,k)  =  b(l  -  b)k_1  k  =  1,2,... 

Note  that  <?b/<?j  =  0  and  s  0  because  of  the  assumed  homogeneity  of  b(j,k). 

The  algorithm  for  planning  the  search  is,  cf.  Eq.  (1), 

p(i)b(l-b)r(l,n,C)  =  max  p( j)b(l-b)r(j ,n_1 

jeJ 

and  if  the  n'th  Inspection  Is  in  field  of  view  £n  =  i  then  search  plan 
5  =  (C'j »^2» •  ♦  •  *£^)  ^  optimal. 

Let's  consider  the  first  few  Inspections  of  three  different  along  orbit 
searches.  In  each  case  we'll  take  J=5  (N=2)  and  p(-2)  +  p(-l)  +  p(0)  +  p ( 1 ) 

+  p ( 2 )  =  0.8  so  p(-2)  =  p(2)  =  0.08,  p(-l)  =  p(l)  =  0.16,  and  p(0)  =  0.32. 

First  let  the  single  glimpse  probability  of  detection  be  high,  b  =  95%. 

Then  the  optimal  search  plan  for  the  first  20  examinations  is  the  sequence 


S*  «  (0,1, -1,2, -2,0, 1,-1, 2, -2, 0,1, -1,2, -2) 


The  probabilities  of  detection  with  this  allocation  of  effort  are  76%, 
79.80%,  and  79.99%  after  the  first  5,  10,  and  15  looks.  Note  that  P[f]  <0.8 
in  this  case.  Since  the  first  five  looks  of  the  optimal  search  plan  are  In 
each  field  of  view  of  the  search  space,  the  existing  search  plans  yield  the 
same  probability  of  detection  after  completion  as  does  the  optimal  one.  This 
should  not  be  at  all  surprising  since  It's  been  assumed  that  you'll  find  the 
satellite  if  you  look  at  it  (l.e.,  b  =  0.95).  The  difference  between  the 
usual  search  plan  and  an  optimal  one  becomes  more  apparent  as  the  probability 
of  detection  decreases  or  the  maximum  amount  of  effort  (K)  increases. 

To  see  this  more  clearly  let  us  reduce  b  to  55%.  The  optimal  search  plan 
for  the  first  21  examinations  is  (both  this  plan  and  the  one  above  appear  to 
be  periodic) 

C*  =  (0,1, -1,0,2, -2,1, -1,0, 2, -2,1, -1,0, 2,-2, 1,-1, 0,2, -2) 

The  cumulative  probabilities  of  detection  following  the  optimal  plan's  alloca¬ 
tion  of  effort  are  47.52%,  65.38%,  73.42%,  and  77.44%  after  5,  10,  15,  and  20 
looks.  If  we  just  look  In  each  cell  of  the  search  space  the  corresponding 
cumulative  probabilities  of  detection  would  be  (for  the  same  total  effort) 

44%,  63.8%,  72.71%,  and  76.72%  respectively.  Once  again,  not  much  difference 
between  the  optimal  plan  and  the  usual  ones.  Of  course  it's  been  assumed  that 
we  have  a  pretty  good  chance  of  detection. 

To  tip  the  balance  towards  the  optimal  plan  consider  the  case  of  a  faint 
object,  say  b  =  0.15.  The  first  twenty-two  stages  In  the  optimal  plan  are 
6*  =  (0,0, 0,0, 0,0, 1,-1, 0,1, -1,0, 1,-1, 0,1, -1,0, 1,-1, 2, -2) 
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Note  the  heavy  concentration  on  that  field  of  view  of  the  search  space  where 
the  a  priori  target  distribution  is  highest.  The  optimal  search  plan  is 
telling  us  that  when  our  chance  of  detection  is  minimal  (b  =  15%)  then  we'd 
best  not  look  where  it's  not  likely  to  be.  The  cumulative  probabilities  of 
detection  following  the  allocation  of  effort  of  the  optimal  plan  are  17.80%, 
28.58%,  36.94%,  and  43.50%  after  5,  10,  15,  and  20  looks.  The  once  in  each 
field  of  view  search  plan  has  only  a  12%  cumulative  probability  of  detec¬ 
tion  upon  completion.  Thus  the  optimal  plan  is  48%  more  efficient  than  the 
usual  one  after  the  customary  expenditure  of  effort.  Repeating  the  usual 
plan  onct  or  twice  yields  cumulative  probabilities  of  detection  upon  com¬ 
pletion  of  22.2%  and  30.87%.  The  optimal  plans  with  equal  expenditure  of 
effort  are  29%  and  20%  more  effective  in  these  cases. 


IV.  ADVANCES 


There  are  two  areas  of  especial  interest  that  we  have  not  dwelled  upon. 
One  concerns  optimal  searches  for  moving  targets.  The  rigorous  theory  that 
supports  optimal  searches  for  fixed  targets  has  not  been  developed  for  this 
case.  Hence,  while  an  experienced  optimal  searcher  may  pursue  his  (or  her) 
intuition  in  such  matters,  certainty  is  lacking.  The  second  domain  concerns 
optimal  searches  for  targets  by  multiple  searcher$--either  moving  or  fixed 
targets  and  colocated  or  separated  searchers.  Space  based  surveillance 
systems  will  yield  yet  another  order  of  complexity  when  non-colocated  and 
moving  searchers  look  for  moving  targets  in  an  asynchronous  fashion.  Progress 
on  these  topics  is  being  made  by  others  and  in  particular,  the  optimal  mul¬ 
tiple  search  for  constant  brightness,  near-stationary  artificial  satellites 
has  been  solved  (reference  5,  see  reference  6  too). 
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